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Abstract —  In  this  work,  we  examine  ground  target  detection  and 
characterization  from  radar  data,  which  incorporates  a  model- 
based  optimization  method  called  dynamic  logic  (DL).  We  apply 
our  methodology  to  a  prototype  airborne  radar  platform  called 
Gotcha,  developed  by  the  Air  Force  Research 
Laboratory/Sensors  Directorate/Automatic  Target  Recognition 
Division  in  recent  years.  The  aircraft  traces  out  a  circular  path 
around  an  area  of  interest,  and  the  onboard,  side-looking  radar 
transmits  and  receives  energy  at  a  constant  pulse  repetition 
frequency,  while  the  main  beam  direction  is  maintained  at  a 
fixed  aim  point  on  the  ground.  Data  collected  during  any 
appropriate  length  arc  of  the  flight  path  can  be  used  to  create 
synthetic  aperture  radar  (SAR)  images  of  the  ground.  The  data 
can  also  be  used  for  ground  moving  target  indication  (GMTI) 
and  provide  Doppler/Range  imagery  of  the  same  ground  area. 
Our  approach  combines  the  computation  of  Range-Doppler 
surfaces  and  a  variable  target  velocity  backprojection  SAR 
method.  Potential  targets  are  detected  using  multiple 
backprojection  images  and  features  are  extracted  using  adaptive 
mixture  models.  We  demonstrate  the  feasibility  of  our  approach 
using  target  truth  information  provided  with  the  Gotcha 
dataset.  We  outline  the  steps  toward  implementing  a 
comprehensive  automatic  target  tracking  solution  based  on 
presented  methodology. 

I.  Introduction 

The  ability  to  conduct  efficient  surveillance  and  tracking  of 
ground  vehicles  and  dismounts  is  essential  for  modern, 
asymmetric  warfare.  This  need  is  emphasized  by  a  number  of 
ongoing  research  programs  devoted  to  developing  systems 
that  exploit  RF  signatures  of  ground  vehicles  and  dismounts. 
Efficient  and  intelligent  automatic  data  processing  algorithms 
are  required  for  real-time  data  analysis,  target 
characterization  and  detection  of  suspicious  events,  such  as 
activities  related  to  lED  emplacement. 

Although  relevant  and  valuable  information  may  be 
collected  by  sensors,  subtle  variations  and  extremely  large 
amounts  of  data  make  it  challenging  to  exploit  that 
information  in  an  efficient  and  effective  manner.  The 
exploitation  challenge  is  compounded  by  the  complex 


operational  environments  where  modern  warfare  occurs,  such 
as  crowded  urban  areas  or  mountainous  terrain. 

To  address  this  urgent  need,  we  are  developing  radar 
exploitation  tools  using  Dynamic  Logic  (DL),  a  general 
framework  for  model-based  signal  processing  and 
optimization  [1].  More  specifically,  we  are  tailoring  DL  to 
the  problem  of  simultaneous  detection,  classification,  and 
tracking  of  ground  targets  with  respect  to  radar  data  collected 
from  a  circular  airborne  path.  Previously,  we  used  DL  to 
develop  algorithms  for  detection  and  tracking  targets  in 
conventional  GMTI  radar  data  having  low  signal-to- 
interference  ratios  with  significant  improvement  [2-6]. 


Figure  1.  Schematic  illustration  of  the  circular  radar  data 
collection. 

In  the  following  sections  we  present  our  general  approach 
to  target  characterization  and  describe  exploitation  using  the 
Gotcha  radar  experimental  data,  which  was  collected  during  a 
series  of  tests  conducted  by  AFRL  at  Wright-Patterson  Air 
Force  base  [7]  [8].  We  then  show  an  example  of  target 
detection  using  DL  and  verify  performance  with  quantitative 
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truth  data.  Finally,  we  discuss  the  results  of  our  experiment 
and  outline  the  directions  of  future  research  in  Gotcha  radar 
characterization  we  intend  to  pursue. 

II.  Radar  Data  Processing 

The  flight  path  of  the  data  collection  aircraft  is  illustrated  in 
Fig.  1.  Side  looking  radar  emits  a  beam  illuminating  an  area 
on  the  ground  with  a  diameter  close  to  1 .5  km,  where  the  slant 
range  to  the  center  of  this  circle  is  about  10  km.  The  radar 
operates  in  spotlight  mode  so  the  same  aim  point  of  the 
antenna  is  maintained  during  the  full  circular  orbit.  The  radar 
characteristics  are  selected  to  achieve  a  region  of 
unambiguous  range  resolution  close  to  1.5  km  long.  The  radar 
uses  stretch  mode  of  operation  over  a  range  of  frequencies  at 
X-band.  Each  returned  pulse  is  stored  in  a  series  of  I  and  Q 
components,  which  constitute  the  raw  data  available  for 
exploitation. 
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Figure  2.  Monitored  area  in  Range/Doppler  space. 

Horizontal  axes  is  Range  in  meters,  vertical  axes  is  Doppler 
in  m/sec. 

SAR  achieves  high  resolution  by  processing  multiple  echo 
waveforms  received  at  different  antenna  positions.  In  our  case 
the  consecutive  antenna  positions  are  located  on  an  arc  of  the 
aircraft  path.  We  select  a  coherent  processing  interval  (CPI) 
containing  Nt  pulses,  each  pulse  containing  Nf  frequencies. 
Thus  the  CPI  data  is  an  Nt  by  Nf  array  of  complex  I  +  Qi 
values.  Standard  Range-Doppler  processing  consists  of 
performing  2-dimensional  Discrete  Fourier  Transform  (DFT) 
of  the  data  for  each  CPI  and  appropriate  scaling  of  the 
units  [10].  An  example  of  such  Range/Doppler  image  from 
the  Gotcha  data  collection  in  2006  is  shown  in  Fig.  2.  The 
range  axis  is  scaled  in  meters  with  the  origin  at  the  antenna 
aimpoint.  The  Doppler  axis  is  scaled  in  meters  per  second  and 
also  has  its  origin  at  the  aimpoint  since  the  radial  velocity  of 
the  aimpoint  relative  to  the  aircraft  is  zero.  All  other  points  on 
the  ground  have  induced  radial  velocity  which  depends  on 
their  location  relative  to  the  aimpoint.  This  induced  velocity 
explains  a  large  central  region  of  the  Range/Doppler  image. 
The  energy  of  the  beam  is  focused  on  the  aimpoint  and 
therefore  weakens  toward  the  edges  of  the  image  showing 


decreased  brightness  of  the  images  of  static  structures 
(buildings,  roads,  etc).  Targets  possessing  radial  velocity 
relative  to  the  aircraft  and  moving  relatively  fast  (at  least  10 
m/sec  in  the  radial  direction)  can  be  detected  since  they  appear 
outside  the  central  clutter  region.  However,the  slowly  moving 
targets  are  located  within  the  central  region.  They  can  be 
detected  by  considering  a  series  of  Range/Doppler  images  and 
identifying  bright  patches  moving  relative  to  the  static 
structures.  Such  a  task  is  complicated  by  the  presence  of  large 
number  of  stationary  objects  in  each  image  and  the  jitter 
induced  by  the  aircraft  motion.  Dynamic  Logic  approach 
described  below  provides  a  good  solution  to  this  problem. 

As  mentioned  above,  jitter  from  air  turbulance  onboard  the 
aircraft  causes  the  antenna  aimpoint  to  randomly  fluctuate 
around  the  desired  location,  which  in  turn  makes  corellations 
between  processed  imagery  and  ground  truth  data  more 
difficult.  A  backprojection  algorithm  [9]  allows  the  formation 
of  more  precise  images  due  to  specifics  of  Gotcha  data 
collection.  This  algorithm  uses  only  the  range  information  in 
the  stored  data.  The  post-processing  done  on  board  of  the  data 
collection  aircraft  includes  motion  compensation  adjusting  the 
phase  of  each  pulse  based  on  the  actual  range  to  the  aimpoint. 
Thus  the  range  information  contained  in  the  data  is  exact. 

This  fact  is  exploited  by  the  backprojection  algorithm  to 
produce  high  resolution  images  of  the  static  structures  on  the 
ground.  Just  like  the  Range/Doppler  imagery,  the  positions  of 
the  moving  targets  are  shifted  due  to  additional  Doppler 
component.  For  example  a  moving  vehicle  will  appear 
moving  in  the  areas  without  roadways.  If  the  radial  velocity 
of  the  vehicle  is  known  the  backprojection  algorithm  can  be 
modified  to  adjust  the  vehicle  position.  The  unknown 
velocity  can  be  estimated  using  the  Range/Doppler  images 
described  above  and  used  to  produce  motion  compensated 
backprojection  images  taking  advantage  of  the  precise 
knowledge  of  the  range.  The  detection  of  the  moving  target 
still  requires  a  good  data  processing  approach  to  separate 
targets  from  static  clutter,  provided  by  Dynamic  Logic 
algorithm  in  our  case. 

As  an  illustration,  consider  a  backprojection  image  in  Fig. 

3.  This  image  is  obtained  by  processing  Nt  =  2000  pulses 
from  Gotcha  data  set  [7].  A  car  stopped  at  the  intersection 
appears  at  its  true  geographical  position.  Fig.  4  shows  a  patch 
of  the  roadway  to  the  left  of  the  intersection  with  the  same  car 
approaching  the  intersection  with  approximately  constant 
speed.  The  car  cannot  be  seen  on  the  image.  We  use  the  GPS 
information  to  estimate  the  velocity  of  the  car  to  adjust  the 
position  of  the  car.  Figure  5  shows  a  sequence  of  4  images 
obtained  by  adjusting  the  backprojection  algorithm  to  the  car’s 
known  velocity.  The  sequence  spans  approximately  1  sec  and 
contains  the  motion  of  the  car  between  the  first  and  the  second 
red  dot  in  Fig.  4  moving  from  left  to  right.  The  level  of  noise 
in  this  image  sequence  underlines  the  need  for  a  robust 
automatic  detection  method. 
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III.  Dynamic  Logic 


Dynamic  Logic  [1]  provides  an  elegant  framework  for 
simultaneously  detecting,  tracking,  and  classifying  multiple 
targets  using  an  adaptive  mixture  of  “competing”  model 
components. 


Figure  3.  Image  of  a  car  stopped  at  an  intersection  obtained 
using  backprojection  algorithm.  The  position  of  the  car  is 
verified  by  GPS  data.  The  axes  are  in  meters  relative  to  the 
antenna  aimpoint. 


Model  parameters  converge  to  the  solution  which 
maximizes  the  similarity  between  the  models  and  the  data, 
providing  track  estimates  and  features  for  target  classification. 


I 


Figure  4.  Backprojection  image  of  the  same  roadway  as 
appears  in  Fig.  3  approaching  the  intersection  from  the  left. 
The  car  is  known  to  move  towards  the  intersection  with 
ground  velocity  ~14  m/sec.  The  red  dots  show  the  GPS 
positions  of  the  car  with  1  sec  increment.  The  car  itself  does 
not  appear  on  this  image. 


The  main  components  of  the  DL  framework  are  the  input  data 
and  the  parametric  models. 


Time  16:48:29+0  sec 


Time  16:48:29+0.23025  sec 


Time  16:48:29+0.4605  sec 


Time  16:48:29+0.921  sec 


Figure  5.  Sequence  of  backprojection  images  of  the  road  in 
Fig.  4  adjusted  by  the  radial  velocity  of  the  vehicle.  The  road 
itself  does  not  appear  in  the  image.  The  rectangle  shows  the 
true  position  of  the  vehicle. 


We  denote  the  input  pixel  by  Xn,  n=l  ..N.  Each  pixel  is  a 
vector  with  components  giving  its  image  coordinates,  color, 
and  timestamp.  There  are  a  total  of  N  data  points  coming 
from  multiple  frames.  We  denote  the  models  by  Mh,  h=l..H, 
where  H  is  the  total  number  of  models.  Each  model  h 
depends  on  parameters  Shi  Mh=Mh(Sh).  Models  Mh  are 
probabilistic,  with  conditional  probability  density  function 
(PDF)  of  the  data  element  Xn  given  the  model  Mh  denoted  by 
l(n|h).  Note  that  this  is  a  function  of  model  parameters.  The 
total  conditional  likelihood  of  the  data  x={  Xn  }  given  all  the 
models  M  =  {  Mh  }  is  expressed  as  follows. 
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Figure  6.  Iterations  of  the  DL  algorithm;  first  column  contains  5  input  frames.  These  images  are  the  same  as  in  Fig.  5 
after  the  application  of  threshold  to  remove  most  of  the  noise.  The  images  are  analyses  by  the  DL  algorithm  to  detect  a 
target  track.  Columns  2,  3,  and  4  illustrate  the  convergence  of  the  algorithm.  The  red  ellipse  is  the  track  model.  The 
clutter  model  is  not  shown.  The  algorithm  converges  to  the  target  track  corresponding  to  GPS  truth  after  10  iterations  . 


L{x  I  M)  = 


n 


h) 


(1) 


Here  rh  is  the  prior  probability  of  encountering  a  model  of 
class  h.  Maximization  of  (1)  with  respect  to  model 
parameters  results  in  the  maximum  likelihood  parameter 
estimate.  This  is  achieved  by  using  the  following  iterative 
procedure.  First,  the  parameters  are  set  to  some  initial  values. 
Second,  the  following  association  weights  are  computed 
using  the  parameter  values 

=  (2) 

S''*''*"  I 

h'=\ 

Next  the  new  parameter  values  are  determined  by  solving  the 
following  equation. 

(3) 

n=\ 


The  algorithm  iterates  between  (2)  and  (3)  until  estimated 
parameters  stop  changing  or  some  other  stopping  criteria  are 
satisfied.  Practical  application  of  this  algorithm  consists  of 
specifying  the  exact  form  of  PDF  (1)  for  each  model  type. 

The  track  model  PDF  consists  of  three  components:  object 
model,  feature  model  and  motion  model.  This  is  expressed  as 
the  following  product  of  PDF’s. 

l{n\h)^  {n  I  h)l {n  \  (n  \  h)  (4) 

Specific  form  of  PDF  depends  on  application.  The  main 
advantage  of  using  this  algorithm  is  that  it  is  more  efficient 
due  to  its  ability  to  avoid  combinatorial  complexity.  This  is 
achieved  by  employing  probabilistic  association  weights  (2) 
and  enforcing  gradual  and  simultaneous  convergence  of 
association  probabilities  and  model  parameters  to  their 
optimal  estimates. 

In  this  work  we  use  the  algorithm  to  detect  moving  targets 
in  sequences  of  images.  Targets  can  be  described  as  blobs 
moving  across  the  image  with  constant  velocity.  The  images 
contain  clutter  in  the  form  of  other  blobs  and  individual 
points  randomly  appearing  is  different  locations.  The 
brightness  of  the  clutter  varies  with  some  blobs  appearing 
brighter  than  the  targets. 
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We  use  linear  motion  model,  Gaussian  feature  model  and  achieve  robust  tracking.  The  combined  algorithm  will  operate 

Gaussian  object  model  given  by  the  following  equations.  as  follows. 


^obj 


V(2^)-|CO. 


^feat 


-0.5(x„-xcf,„f  (COf,)  \x„-xcf,„) 


yl{2^y\CX, 


-0.5{xc,„-X,„  f  (CX, )-'  (xc,„-X,„ ) 


K5) 


The  motion  is  incorporated  in  the  motion  PDF  using  the 
following  equation. 

X,„^X0,+V,t„  (6) 

Here  the  model  parameters  are  the  mean  and  covariance  of 
each  of  the  models  and  the  initial  position  and  velocity  of  the 
track.  The  target  PDF  is  therefore  obtained  by  substituting  (5) 
into  (4).  The  clutter  model  uses  Gaussian  feature  model. 
There  are  no  object  and  motion  models  for  clutter,  so  we 
model  the  position  of  the  clutter  points  with  uniform  density. 
Thus  the  clutter  model  PDF  is  given  as  follows. 


l{n  I  c)  = 


i _ ^  ^-O-HYn-Ycf  {CY,)-\y,-YJ 

^  p)rY\cy^\ 


(7) 


Here  A  is  the  area  of  the  image.  The  dimensionality  of  the 
feature  d=l  in  our  case  and  we  use  model  index  c  instead  of  h 
to  emphasize  that  this  is  the  clutter  model  PDF. 

We  applied  the  algorithm  to  the  sequence  of  images  in  Fig 
.6.  The  images  were  thresholded  to  get  rid  of  most  but  not  all 
background  noise.  Thresholding  is  not  necessary  for 
application  of  the  algorithm  and  was  done  to  speed  up  the 
processing.  The  algorithm  quickly  converged  to  the  target 
track  separating  image  points  coming  from  the  target  from 
clutter  points.  The  run  time  of  this  example,  implemented  in 
Matlab,  is  2.7  seconds  on  a  Dell  desktop  computer  with 
3GHz  processor  and  4GB  of  RAM.  This  can  be  improved  by 
several  orders  of  magnitude  with  more  efficient 
implementation  and  specialized  hardware. 


IV.  Discussion 

The  goal  of  this  work  was  to  describe  our  current  approach  to 
the  Gotcha  data  exploitation.  We  discussed  two  ways  of 
forming  SAR  imagery  from  the  data.  We  gave  a  brief 
overview  and  an  example  of  usage  of  the  Dynamic  Logic 
algorithm  for  detection  and  feature  extraction  from  the  SAR 
imagery.  The  long  term  purpose  of  this  research  is  the 
development  of  comprehensive  solution  for  detection  and 
tracking  of  moving  targets  -  vehicles  and  dismounts.  We 
propose  to  combine  the  Range/Doppler  and  the  backprojection 
SAR  imagery  to  exploit  the  advantages  of  both  methods  to 


1)  The  ground  radial  speed  of  the  moving  target  is  estimated 
using  DL  algorithm  applied  to  Range/Doppler  images.  The 
Range/Doppler  images  can  be  computed  relatively  fast. 

2)  The  estimated  target  speed  is  plugged  into  the 
backprojection  algorithm  to  obtain  the  SAR  images  of  the 
moving  target.  The  precision  of  backprojection  images  results 
in  exact  positioning  of  the  target. 

3)  Use  DL  algorithm  to  detect  and  estimate  the  parameters  of 
the  track  and  other  features  of  the  target  such  as  its  radar  cross 
section. 

4)  Once  the  target  is  detected  the  subsequent  tracking  can  be 
done  by  extending  the  track  over  the  future  backprojection 
images  and  using  DL  and  parameters  estimated  during 
detection. 
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